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Abstract 

We examine the thermodynamic stability of large black holes in four-dimensional anti- 
de Sitter space, and we demonstrate numerically that black holes which lack local 
thermodynamic stability often also lack stability against small perturbations. This 
shows that no-hair theorems do not apply in anti-de Sitter space. A heuristic argument, 
based on thermodynamics only, suggests that if there are any violations of Cosmic 
Censorship in the evolution of unstable black holes in anti-de Sitter space, they are 
beyond the reach of a perturbative analysis. 
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1 Introduction 



The Gregory-Laflamme instability [1] is a classical instability of black brane solutions 
in which the mass tends to clump together non-uniformly. The intuitive explanation 
for this instability is that the entropy of an array of black holes is higher for a given 
mass than the entropy of the uniform black brane. The intuitive explanation leaves 
something to be desired, since it applies equally to near-extremal Dp-branes: scaling 
arguments establish that a sparse array of large black holes threaded by an extremal 
Dp-brane will be entropically favored over a uniform non-extremal Dp-brane; however it 
is not expected that near-extremal Dp-branes exhibit the type of instability found in [1]. 
It was checked in [2] that a Dp-brane which is far from extremality (that is, one whose 
tension is many times the extremal tension) does have an instability. It was also shown 
that the instability persists for charged black strings in five dimensions fairly close to 
extremality. 1 Less is known about the case of near-extremal D3-branes, M2-branes, and 
M5-branes, but one may take the absence of tachyons in the extensive AdS-glueball 
calculations ([3, 4] and subsequent works — see [5] for a review) as provisional evidence 
that these near-extremal branes are (locally) stable. 2 

In the AdS/CFT correspondence [6, 7, 8], one might at first think that the existence 
of a unitary field theory dual forbids an instability. But suppose we are at finite 
temperature, and that there is a thermodynamic instability in the field theory — like 
the onset of a phase transition. Then it is quite natural for some fluctuation mode (or 
modes) to grow exponentially in time, at least in a linearized analysis, as one nucleates 
the new phase. Exciting an unstable mode is a change in the state of the field theory, 
not its lagrangian; thus according to AdS/CFT there should be a normalizable mode 
in AdS which likewise grows exponentially with time [9]. This might be referred to as a 
"boundary tachyon," or a "tachyonic glueball," since in the gauge theory it corresponds 
to some bound state with negative mass-squared. We will prefer the term "dynamical 
instability," which is meant to convey that there is an instability in the Lorentzian time 
evolution of the black brane, in both its supergravity and dual field theory descriptions. 

To sum up, the existence of a field theory dual makes plausible the following adap- 

1 The charged black string studied in [2] happens to be thermodynamically unstable all the way 
down to extremality: the specific heat is negative. Thus (1) would lead us to believe that this non- 
extremal black string is always unstable. The extremal solution should be stable since it can be 
embedded in a supersymmetric theory as a BPS object. 

2 More properly we should say that the near-extremal black brane solutions with many units of D3- 
brane, M2-brane, or M5-brane charge appear to be stable. A single brane has Planck scale curvatures 
near the horizon, so classical two-derivative gravity does not provide a reliable description. We will 
concern ourselves exclusively with solutions which have a discrete parameter (M2-brane charge, for 
the most part) which can be dialed to infinity to suppress all corrections to classical gravity. 
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tation of the entropic justification for the Gregory-Laflamme instability: 

For a black brane solution to be free of dynamical instabilities, it is necessary 
and sufficient for it to be locally thermodynamically stable. 

Here, local thermodynamic stability is defined as having an entropy which is concave 
down as a function of the mass and the conserved charges. This criterion was first used 
in a black brane context in [10], where it was found that spinning D3-branes could 
be made locally thermodynamically unstable if the ratio of the spin to the entropy 
was high enough. Further work in this direction, relevant to the current paper, has 
appeared in [11, 12, 13, 14]. For a somewhat complementary point of view on the 
nature of the unstable solutions, see [15, 16]. 

The conjecture (1) is meant to be a local version of the argument about whether 
the array of black holes or the black brane has higher entropy; however it seems on 
more precarious ground since one may not be able to write down a non-uniform sta- 
tionary solution that competes with the black brane entropically. Nonetheless, it was 
shown in [17] that (1) predicts with good accuracy the value of the charge where the 
four-dimensional anti-de Sitter Reissner-Nordstrom solution (AdS^-RN) develops an 
instability. 

The aim of this paper is to give a fuller exposition of the calculations in [17], to 
present a more complete picture of the results of the numerics, and to explore via ther- 
modynamic arguments the likely paths for time-evolution of the unstable solutions. 
Section 2 contains a summary of the AdS^-KN solution and some generalizations of 
it in Af = 8 gauged supergravity and in higher dimensions. Section 3 discusses the 
thermodynamic instability which occurs for large charge. In section 4, a linear pertur- 
bation analysis is carried out around the AdS^TW solution. In the large black hole 
limit, a dynamical instability appears when local thermodynamic stability is lost. The 
existence of a dynamical instability was the main result of [17]. It disproves the claim 
of [18, 19] that charged black holes in AdS are classically stable. As we explain in 
section 4, the instability persists some ways away from the large black hole limit, pro- 
viding the first proven example of a black hole with a compact horizon and a pointlike 
singularity which exhibits a dynamical Gregory-Laflamme instability. 3 Such solutions 
are interesting from the point of view of Cosmic Censorship, and we discuss the possi- 
bility of forming a naked singularity, or at least regions of arbitrarily large curvatures. 
Our main result here is that adiabatic evolution toward maximum entropy does not 
lead to solutions which arise from making the mass smaller than some appropriate 

3 Here we are referring to the existence of a local instability visible in a classical analysis. It has 
been observed [20] that the AdS-Schwarzschild solution times a sphere can have a lower entropy than 
a Schwarzschild black hole of the same mass which is localized on the sphere. This demonstrates 
global but not local instability, and suggests the possibility of tunneling from one configuration to the 
other. 
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combination of the charges. Because entropic arguments appear to give good informa- 
tion not only on the existence of dynamical instabilities but also on the direction they 
point, it is reasonable to predict from our results that no perturbative analysis of a 
smooth black hole in AdS will demonstrate a violation of Cosmic Censorship. 

The unstable mode of the AdS^-KN solution does not involve fluctuations of metric 
at linear order. Rather, it involves the gauge fields and scalars of Af = 8 gauged 
supergravity. Because the metric is not fluctuating, it may seem odd to describe the 
process as a Gregory-Laflamme instability. But we claim that the instability we see 
is in the same "universality class" as instabilities where the horizon does fluctuate: 
to be more precise, if the charges of the black hole are made slightly unequal, then 
generically the instability will involve the metric. In fact, the metric does fluctuate in 
the equal charge case as well — only at a subleading order that is beyond the scope of our 
linearized perturbation analysis. We would in fact make the case that any dynamical 
instability of a black hole which leads to non-uniformities in charge or mass densities 
should be considered in the same category as the Gregory-Laflamme instability of 
uncharged black branes. 

We emphasize that this paper is concerned with the relation between local thermo- 
dynamic stability of stationary solutions and the stability of their classical evolution 
in Lorentzian time. It is known [21, 22, 23, 24, 25, 26] that black holes which are 
thermodynamically unstable have an unstable mode in the Euclidean time formalism. 
For spherically symmetric black holes this mode is an s-wave. The interpretation is 
that, for instance, an AdS-Schwarzschild black hole in contact with a thermal bath 
of radiation will not equilibrate with the bath if the specific heat of the black hole is 
negative. This beautiful story does not fall under the rubric of problems we are con- 
sidering, because the processes by which equilibration takes place in Lorentzian time 
include Hawking radiation, which is non-classical. Rather, we are contemplating black 
holes or branes in isolation from other matter, in a classical limit where Hawking radi- 
ation is suppressed, and inquiring whether a stationary, uniform black object wants to 
stay uniform or get lumpy as Lorentzian time passes. It is less clear that there should 
be any relation between this dynamical question and local thermodynamic stability: 
for instance, a Schwarzschild black hole in asymptotically flat space is stable. 4 Yet we 
conjecture that (1) gives a precise relation when the black object has a non-compact 
translational symmetry. 

Our focus in this paper is black holes in AdS and their black brane limits; however 
the conjecture (1) is intended to apply equally to any black brane. It may apply 
even beyond the regime of validity of classical gravity. Any "sensible" gravitational 

4 This stability is implied by classical no-hair theorems, see for example [27]. A more extensive list 
of references on no-hair theorems can be found in [28]. A consequence of the present work is that 
these theorems cannot be extended to charged black holes in AdS. 
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dynamics should satisfy the Second Law of Thermodynamics, and (1) is motivated 
solely by intuition that Lorentzian time evolution should proceed so as to increase 
the entropy. (The stipulation of translational invariance prevents finite volume effects 
from vitiating simple thermodynamic arguments). For instance, it has recently been 
shown [29] that the near-extremal NS5-brane has a negative specific heat arising from 
genus one contributions on the string worldsheet (see also [30, 31], and [32] for related 
phenomena in 1+1-dimensional string theory). 5 This is not classical gravity, but (1) 
leads us to expect an instability in the Lorentzian time evolution of near-extremal 
NS5-branes. 6 The instability would drive the NS5-brane to a state in which the energy 
density is non-uniformly distributed over the world-volume. 

2 The Ad^-RN solution and its cousins 

The bosonic part of the lagrangian for M = 8 gauged supergravity [33, 34] in four 
dimensions involves the graviton, 28 gauge bosons in the adjoint of 5*0(8), and 70 
real scalars. Because of the scalar potential introduced by the gauging procedure, flat 
Minkowski space is not a vacuum solution of the theory; rather, AdS^ is. It is known 
[35] that the maximally supersymmetric AdS± vacuum of Af = 8 gauged supergravity 
represents a consistent truncation of 11-dimensional supergravity compactified on S 7 . 
The AdS^ x S 7 solution can be obtained as the analytic completion of the near-horizon 
limit of a large number of coincident M2-branes. 7 Making the M2-branes near-extremal 
corresponds to changing AdS^ to the Aci^-Schwarzschild solution. Near-extremal M2- 
branes can also be given angular momentum in the eight transverse dimensions. There 
are four independent angular momenta, corresponding to the U{1) A Cartan subgroup 
of SO (8): these reduce to electric charges in the AdS 4 description. The electrically 
charged black hole solutions can be obtained most efficiently by first making a consis- 
tent truncation of the full Af = 8 gauged supergravity theory to the U{1) A gauge fields 
plus three real scalars. Consistent truncation means that any solution of the reduced 
theory can be embedded in the full theory, with no approximations. For our purposes, 
it can be viewed as a sophisticated technique for generating solutions. The truncated 

5 We thank D. Kutasov for bringing [29, 32] to our attention. 

6 We thank M. Rangamani for a number of discussions on this point. 

7 As stated in the introduction, taking the number of M2-branes large makes the geometry smooth 
on the Planck/string scale and thus suppresses corrections to classical two-derivative gravity. 
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where the signs on the gauge fields can be chosen independently. We will lose nothing 
by choosing them all to be +. The quantities Qa are the physical conserved charges, 
and they correspond to the four independent angular momenta of M2-branes in eleven 
dimensions. The mass is [11] 
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where zh is the largest root of F(zh) = 0. Only for sufficiently large At do roots to this 
equation exist at all. When they don't, the solution is nakedly singular. 

We will be most interested in the case where all four charges are equal, qA = q. 
Then the solution can be written more conveniently in terms of a new radial variable, 
r = z + q, and it takes the form 
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with the scalars set to 0. In (6), F 0r is the common value of all four gauge field 
strengths F^\ The geometry (6) is a solution of pure Einstein-Maxwell theory with 
a cosmological constant: it is the AdS^-KN solution. 

There are related solutions to maximally supersymmetric gauged supergravity in five 
and seven dimensions, corresponding respectively to spinning D3-branes and spinning 
M5-branes. In the case of D3-branes, there are six transverse dimensions, the rota- 
tion group is 5*0(6), the Cartan subalgebra is U(l) 3 , and as a result there are three 
independent angular momenta (or charges in the Kaluza-Klein reduced description). 
In the case of M5-branes, there are five transverse dimensions, the rotation group is 
SO (5), the Cartan subalgebra is U(l) 2 , and there are two independent angular mo- 
menta/charges. We will record here only the Einstein frame metric in the Kaluza-Klein 
reduced description, in conventions where Gn = 1 and L is the radius of the asymp- 
totic AdS space. For further information on these solutions, the reader is referred to 
[11, 15]. The metrics are 

AdS 5 : ds 2 = -H^Fdt 2 + H*(^- + r 2 d^ 

H=f[h A F = l-± + £ H h A = l + qA 



„2 

A=l 



L 2 ' r- 



AdS 7 : ds 2 = -H-lFdt 2 + H l *(^ + r 2 dn^ (7) 

2 2 

H=l[hA F = l-^+ r -H h A = l + f A . 

A=l r 



3 Thermodynamics 
3.1 Generalities 

Given the solutions (3) and (7), we may read off the entropy, the mass, and the con- 
served electric charges. Typically it is most straightforward to express these quantities 
in terms of the non-extremality parameter ji and the boost parameters Pa- However it 
is possible to eliminate \i and f3 A and find a polynomial equation relating M, S, and 
the Q A . This equation can be solved straightfowardly for M, but not in general for 
S. We will quote explicit results for M = M(S,Qi, . . . ,Q n ) in the next section. In 
this section we will discuss thermodynamic stability assuming that M(S, Q±, . . . , Q n ) 
is known. 

The micro canonical ensemble is usually specified by a function S = S(M, Qi, . . . , Qn)- 
Assuming positive temperature (which is safe for regular black holes since the Hawk- 
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ing temperature can never be negative), one may always invert M = M(S,Qa) to 
S = S(M, Qa), where now we abbreviate Qi, . . . , Q n to Qa- A standard claim in clas- 
sical thermodynamics is that the entropy for "sensible" matter must be concave down 
as a function of the other extensive variables. Locally this means that the Hessian 
matrix, 

/ a 2 s d 2 S \ 

ttS - dM 1 OMOQb \ fo\ 

m,q a = I d 2 s a 2 s I ' y°) 

\ 8Qa9M 8Qa9Qb > 

satisfies Hf f g A < 0, i.e. it has no positive eigenvalues. To understand what this 
requirement means, consider the simplest case where n = and d 2 S/dM 2 > 0. This 
is the statement that the specific heat is negative. A substance with this property (in 
a non-gravitational setting, but equating mass with energy) is unstable: if we start 
at temperature T, then it is possible to raise the entropy without changing the total 
energy by having some regions at temperature T + ST and others at T — ST. Since we 
are implicitly assuming a thermodynamic limit, it is irrelevant how big the domains of 
high and low temperature are. In a more refined description (e.g. Landau- Ginzburg 
theory), these domains might have a preferred size, or at least a minimal size. 

In the more general setting of many independent thermodynamic variables, let us 
define intensive quantities 

(2/o, y u ..., y n ) = (M/V, Q 1 /V, Q n /V) , (9) 

where V is the volume. Suppose that Hf^g^ has a positive eigenvector: ~Rm q a v = \v 
with A > 0. Through a variation 

Vi -»■ Vj + ev i > ( 10 ) 

where e is a function of position which integrates to 0, we can raise the entropy without 
changing the total energy or the conserved charges. Thus positive eigenvectors of 
H^q^ indicate the way in which mass density and charge density tend to clump. 
Presumably the eigenvector with the most positive eigenvalue gives the dominant effect. 

The stability requirement Hf 1 Q A < may be rephrased as H^q a > 0, where H^q a 
is the Hessian of M with respect to S and Qa- This is easy to understand from a 
geometrical point of view. Hm q a < says that all the principle curvatures of S(M, Qa) 
point toward negative S, or, equivalently, away from the point (S, M, Qa) = (oo, 0, 0). 
Now, the point (S, M, Q A ) = (0, oo, 0) is on the opposite side of the co-dimension 
hypersurface defined by S = S(M, Qa) from (S, M, Qa) = (0, oo, 0). Thus all principle 
curvatures should point toward (0, oo,0), which means that Hf t Q A > 0. To determine 
the region of thermodynamic stability we may thus require det H^q a > 0, and then 
take the smallest connected components around points which are known to be stable. 

While regions of stability are conveniently calculated from H^q a , it is not clear that 
the eigenvector of Hm q a with the largest positive eigenvalue can be read off easily from 
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Hg t Q A . So it is useful to express "H.mq a directly in terms of derivatives of M(S, Qa)'- 
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A prescription for dealing with energy functions which violate the convexity condition 
Hg t Q A < is the Maxwell construction, where one replaces M(S, Qa) with its convex 
hull (or S(M,Q A ) by its convex hull — it's the same thing). This formal procedure is 
equivalent to allowing mixed phases where some domains have higher mass density or 
charge density than others. The energy functions resulting from charged black holes in 
AdS have the curious property that the convex hull is completely flat in some directions, 
so that chemical potentials (after taking the convex hull) are everywhere zero. This 
arises because, in certain directions, M rises slower than any nontrivial linear function 
of the other extensive variables. In this situation the Maxwell construction does not 
make much sense, because the mixed phases that it calls for have charges and mass 
concentrated arbitrarily highly in a small region, while the rest of the "sample" is 
at very low charge and mass density. A similar example in the simpler context of 
no conserved charges would be a mass function M(S) like the one in figure 1. Here 
the natural physical interpretation is that the region between A and B represents a 
stable phase, while the region to the right of B is unstable toward clumping most 
of its energy into small regions. This tendency would presumably be cut off by some 
minimal length scale of domains. The mass functions obtained from charged black holes 
in AdS look roughly like figure 1 along some slices of the space of possible (S, Qa)- The 
interpretation we will offer is that the black holes are stable in the regime of parameters 
where convexity holds, and that they become dynamically unstable toward clumping 
their charge and energy outside this region. 8 

The line of thought summarized in the previous paragraph was already advanced in 
[12], but with only thermodynamic arguments to support it. A competing point of view 
was suggested in [18]: the black holes in question have no ergosphere (more precisely, 
there is Killing vector field which is timelike everywhere outside the horizon), and this 

8 There is a subtlety, discussed in [12], about the precise location of the boundary between stable 
and unstable regions. As the system approaches the inflection point at B, finite fluctuations might 
allow it to make small excursions into the unstable region. Working in a large N limit where classical 
supergravity applies on the AdS side of the duality seems to suppress such fluctuations. 
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was argued to imply that there could be no superradiant modes, and hence no classical 
instability in the Lorentzian-time dynamics. The argument used the dominant energy 
condition, which need not always be satisfied by matter in AdS: in fact, the scalars 
tfi in (2) violate the dominant energy condition because of their tachyonic potential 
(which however does satisfy the Breitenlohner-Freedman bound). 

In [17], an explicit numerical calculation demonstrated the existence of a dynamical 
instability for certain AdS^-TW black holes (related to spinning M2-branes with all four 
spins equal, as explained in the previous section). We will discuss this calculation at 
greater length in section 4. For now let us only remark that in the limit of large black 
holes, where the horizon area is infinite, the instability appears when thermodynamic 
stability is lost, up to a discrepancy of 0.7% which we suspect is numerical error. 
Furthermore, the combination of supergravity fields which became unstable indicated 
a change in local charge densities precisely in agreement with the analysis leading 
up to (10). Thus the conjecture (1) was tested to reasonably good accuracy along a 
two-parameter subspace (entropy and the common value of the four charges) of the five- 
parameter phase space. Further tests in AdS^ are significantly more difficult because 
the metric usually enters in to the perturbation equations in a non-trivial way. However 
we will indicate in section 4 another case where the metric decouples. Tests in AdS§ 
and AdSj can also be performed most easily in the equal charge case, but the analysis 
is somewhat more tedious because the spinor formalism is not as well worked out in 
higher dimensions (and probably is more cumbersome in any case). 

Despite the absence of comprehensive tests, we will use (1) and the idea that black 
hole perturbations should follow the most unstable eigenvector of Hmq a to propose 
in section 3.3 a qualitative picture of the evolution of unstable black holes in AdS. In 
brief, once the boundary of stability is passed, the independent charges tend to clump 
separately, as if they repelled one another but attracted themselves. But this is only 
an approximate tendency, with significant exceptions to be noted in section 3.3. When 
a particular angular momentum density becomes very large, it is possible that anti-de 
Sitter space fragments by a classical process. 9 This also will be discussed at greater 
length in section 3.3. We emphasize that our proposals for the evolution of unstable 
black holes are largely conjectural, and difficult to check by any means other than 
numerical solution of the full equations of motion. 

fragmentation of AdS via tunneling has been discussed in [37]. 
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3.2 Explicit formulas 

It is possible to eliminate all the auxiliary quantities from (3), (4), and (5), and express 
M directly in terms of the entropy and the physical charges as 

i 

(12) 

We will often be interested in the limit of large black holes, M/L ^> 1. In this limit 
we have 

i 

(13) 

with corrections suppressed by powers of M/L. As M/L approaches infinity, one 
obtains a black brane solution in the Poincare patch of AdS±. Formally this limit can 
be taken by expanding (3) to leading order in small fa, dropping the 1 from F, and 
replacing S 2 by R 2 in the metric. 

As remarked in the previous section, local thermodynamic instability can be ex- 
pressed as convexity of the function M{S,Q\,Q2,Qz,Qi)- By setting the Hessian of 
(13) equal to zero, we obtain the boundary separating the stable from the unstable 
region: 

3S 8 - 2n 2 L 2 S 6 Y,Q\ + tt 4 L 4 S 4 £ (Q a Qb) 2 - vr 8 L 8 {[Q 2 A = 0. (14) 

A=l A<B A=l 

First let us consider the case where the charges are pairwise set equal: Qi = Q3 and 
Q2 = Qa- The above equation then factorises, giving us three relevant factors: 

(S 2 - n 2 L 2 Ql) (S 2 - n 2 L 2 Q 2 2 ) -^f(Ql + Ql + y/Qf + Qt + UQ 2 Q 2 )j = . 

(15) 

When at least one of these factors become negative, H^q a develops a negative eigen- 
vector and the black hole becomes thermodynamically unstable. A more convenient 
form may be obtained by eliminating S in favor of the mass M and introducing the 
dimensionless variable Ys = — rH-- The above three equations in the new variables 
become 

xi + xl + Sxlxi + (xl + xl) ^x\ + xi + i±xlxi\ 2 

- 54 (x\ + xl + \M + Xi + 14x? Xl) = (16) 
X? + 2x?X2(x? + X2)-4 = 

x 6 2 + 2xixi(xl + xi)-4 = o. 
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The region depicting thermodynamically stable black holes is the intersection of the 
areas under the 3 curves as shown in figure 2(a). 

The other relevant curve is the one separating nakedly singular solutions from regular 
black holes. The mathematical criterion for having a regular black hole solution is that 
the polynomial F in (3) should have a zero. In the large black hole limit, and in terms 
of xi an d X2, this criterion reduces to 

xlxKxt + xl) - IxlxKxt + Xi) + l^xlxKxl + Xi) - 4(x? + Xl) + Qxtxt - 432 = . 

_ ^ (17) 

To determine if a black hole with given values of mass and charges is unstable, one 
first computes the values of xi an d X2 and locates this point in figure 2(a). The black 
hole is unstable if the point lies outside the shaded region depicting stable black holes 
but is within the boundary which separates black holes with naked singularities from 
those with a horizon. If the point lies in the unshaded (unstable) region of the plot 
without the vector field shown, it means that within each pair one charge wants to 
increase while the other decreases. The unstable eigenvector has no components along 
the hyperplane Qi = Q3 and Q2 = Qa and is not shown. 

Finally, let us collect the thermodynamic results for the special case of all charges 
equal. We see that thermodynamic instability is present in the narrow region 1 < 
X < v / 3/2 2 / 3 . The associated eigenvector has the form (0,1,-1,1,-1) where the 
components are along the axes M, Q 1 , Q 2 , Q3, and Q4 respectively: it looks like one 
pair of charges wants to increase while the other decreases. This can happen only 
locally, with each of the four charges conserved globally. 

We'll also consider the case in which only two of the charges, Q\ and Q2, are non- 
zero. To get the region of thermodynamically stable black holes, we set Q 3 = Q4 = 
in (14): 

3S 4 - 2ir 2 L 2 S 2 (Qj + Q\) + ix A L^Q\Q\ = . (18) 

Just as we did in the previous case, we first eliminate S in favor of the mass M and 
then introduce the dimensionless variables Xi = ?' 1 to get: 

nxl + xl) + 21xixl(xl + xl) 

(19) 

+ (lOxt + 10x2 + Wxlxl) Vxf + xi - xlxl - 432 = . 

This is the boundary of the stable region, and is plotted in figure 2(b). Unlike the case 
of charges set equal pair-wise, black holes with two charges set to zero always have a 
horizon. This may be connected with the fact that there is a limit of rotating M2-branes 
with only two independent angular momenta nonzero which is a well-defined multi- 
center M2-brane solution, while with all angular momenta nonzero the corresponding 
limit is a singular configuration in eleven dimensions [38, 39, 40]. 
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For black holes in AdS 5 and AdS 7 , we will simply record here the mass in terms of 
the entropy and charges: 



AdS 5 : M 
AdS 7 : M 



2L 2 (2tt 4 5) 2 / 3 
5 

4L 2 (4tt 9 ^) 2 / 5 



nV3 



II (4^ 2 + n 4 L 2 Q 2 A 

4=1 

n (IQS 2 + n 6 L 2 Q 2 A ) 



.A=l 



2/5 



(20) 



Stability analyses similar to the AdS 4 case can be carried out for AdS 5 and AdS 7 . 
Some work along these lines was presented in [12], but the explicit expressions in (20) 
make the calculations much easier. 



3.3 Adiabatic evolution 

Tracking the evolution of unstable black holes in Lorentzian time is difficult. We have 
succeeded in establishing perturbatively the existence of a dynamical instability for 
the very special case of all charges equal: this is explained in section 4. This simplest 
case required the numerical solution of a fourth order ordinary differential equation 
with constraints at the horizon of the black hole and the boundary of AdS&. Most 
other cases for black holes in AdS± involve fluctuations of the metric, which makes 
the analysis significantly harder. To investigate the instabilities beyond perturbation 
theory would require extensive numerical investigation of the second order PDE's that 
comprise the equations of motion of Af = 8 gauged supergravity. 

The aim of this section is to use thermodynamic arguments to guess the qualitative 
features of the evolution of unstable black holes. Here we focus exclusively on the large 
black hole limit; however the conclusions may remain valid to an extent for finite size 
black holes with dynamical instabilities. The intuition is that knowing the entropy as 
a function of the other extensive parameters amounts to knowing the zero-derivative 
terms in an effective Landau- Ginzburg theory of the black hole (or of its dual field 
theory representation). 

As explained in the paragraph around (9) and (10), an unstable eigenvector of H^q^ 
(by which we mean one with positive eigenvalue) suggests a direction in which a black 
hole solution can be perturbed in order to raise entropy while keeping its total mass 
and conserved charges fixed; moreover it was shown in [17] (as we will explain in 
section 4) that the black hole's dynamical instability causes it to evolve in precisely 
the direction that the eigenvector indicates. The physics has no infrared cutoff, as is 
typical in Gregory-Laflamme setups, so we may hope that the charge and mass densities 
vary over long enough distance scales that we may continue to use the most unstable 
eigenvector of Hm ; q a locally to determine the direction of the subsequent evolution. 
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Following this line of thought to its logical conclusion leads us to the claim that the mass 
density and charge densities will locally evolve, subject to the constraints of conserving 
total energy and charge, from their initial values to values along a characteristic curve 
of the unstable vector field of Hf I n A . This can only be approximately correct: finite 
wavelength distortions will occur, and it is not precisely right anyway to say that the 
time-evolution of Einstein's equations proceeds so as to maximize black hole entropy. 
Nevertheless it seems to us likely that a correct qualitative picture will emerge from 
tracking the flows generated by the most unstable eigenvector of Hfj q a . At late times, 
or when charge and mass density are highly concentrated in small regions, another 
description is needed. 

The characteristic curves of the most unstable eigenvector of Hq Ma may terminate 
in a region of stability, or in a region of naked singularities. Cosmic Censorship plus the 
conjectures of the previous paragraph suggest that the latter should never happen. This 
can be checked explicitly for the examples that we have. To this end, one can choose a 
generic value of charges and mass so that the black hole is almost naked, then determine 
the most unstable eigenvector of Hf f g A , and then check that it is tangent to the surface 
separating naked singularities from regular black holes. We carried this out numerically 
for several cases and verified tangency; however we do not have a general argument. It 
appears, in fact, that the normal vector to the surface separating naked singularities 
from black holes is a stable eigenvector of Hf I n A (i.e. its eigenvalue is negative) — at 
least in the three-dimensional subspace with Qi = Q 3 and Qi = Q4 — so the obvious 
approach to an analytic demonstration that Cosmic Censorship is not violated by 
adiabatic evolution of black holes is to show that this normal vector is always a stable 
eigenvector of Hf I g A . For now we content ourselves with the observation that in all 
the cases we have checked numerically, adiabatic evolution does stay in the region of 
regular black holes. 

It is also possible that a characteristic curve becomes unstable at some point, in the 
sense that nearby characteristic curves diverge from it. To refine our previous claim, 
we may suppose that the black hole evolves along a bundle of nearby characteristic 
curves emanating from the original mass and charge density. This bundle may remain 
nearly one- dimensional, or it may split or become higher dimensional. We will not 
investigate the stability properties of the characteristic curves in any detail. Note that 
we are not attempting to specify any spatial or temporal properties of the evolution, 
only the range of mass and charge densities which form. 

We present in figure 2 plots of unstable eigenvectors of the Hessian matrix HfjQ A , 
projected onto a plane parametrized by two of the charges. From these vector fields, 
we may conclude that the different charges exhibit some tendency to separate from one 
another, but that this does not always happen, as in the upper right part of figure 2(b). 
The crucial point is that the unstable eigenvectors don't have a component normal 
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to the boundary between naked singularities and regular black holes. Although this 
appears obvious from figure 2(a), the plot is slightly misleading in that the eigenvectors 
have been projected onto the plane of Q\ = Q3 and Q2 = Qa- One must preserve the 
components of vectors in the M direction to verify tangency. 

When some angular momenta become large compared to the entropy for a spinning 
M2-brane solution, the geometry in eleven dimensions is approximately given by a 
rotating multi-center brane solution [38] . If one angular momentum is large, this multi- 
center solution is in the shape of a disk; if two are large and equal, it has the shape 
of a filled three-sphere. It seems clear that solutions of this form in an asymptotically 
flat eleven-dimensional spacetime are unstable toward fragmentation in the directions 
transverse to the M2-brane. This would mean that anti-de Sitter space fragments. In 
terms of the SU(N) gauge theory, the disk corresponds to a [/(l)^ 1 Higgsing, and 
in the fragmentation process some groups of f/(l)'s try to come together to partially 
restore gauge invariance. It is not certain that such fragmentation occurs, particularly 
if the angular momentum density is large only locally. We merely indicate it as a 
possibility in the complicated late-time evolution of unstable black holes. 

Finally, it is possible that the unstable black holes evolve to a some new stationary 
solution, presumably with a non-uniform horizon. No such solutions are known. It 
seems to us most likely that, upon becoming unstable, black holes in AdS undergo an 
evolution which eventually produces large curvatures. 

4 Existence of a dynamical instability 

The existence of dynamical instabilities for thermodynamically unstable black branes 
should be completely generic. However, as mentioned already, the stability analysis is 
technically complicated for the general case of unequal charges: perturbations of the 
metric, four gauge fields, and three scalars lead to difficult coupled partial differential 
equations. Here we focus on the AdS^-KN example, where the metric decouples and the 
problem can be reduced to a single gauge field and a single scalar. A formal argument 
relating thermodynamic and dynamical instability was suggested in [17], using the 
identification of the free energy with the Euclidean supergravity action; however we 
have not yet succeeded in making this argument rigorous. 

Because the unstable eigenvectors of Hmq a (for all charges equal and sufficiently 
large) do not involve any change in the mass density, it is natural to expect that the 
perturbations that give rise to an unstable mode do not involve the metric. 10 More 
precisely, because of the form of the unstable eigenvectors, we expect that a relevant 

10 Indeed, we suspect that the decoupling of the metric is possible precisely when there is an eigen- 
value of q a which does not have a component in the M direction. 
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perturbation is 

SF A = a\5F (21) 

for some SF and fixed i, where the a l A were defined in (2). In section 3 we saw explicitly 
that 8Qi = 5Qs = —SQ2 = —SQ4 gave an unstable eigenvector; now we make a trivial 
alteration and focus on SQ\ = SQ2 = —SQ3 = —SQ4. Correspondingly we set i — 1 in 
(21). 

The spectrum of linear perturbations to charged black holes in AdS has been consid- 
ered before [41], but for the most part the perturbations under study were minimally 
coupled scalars. It is impractical to sift through the entire spectrum of supergravity 
looking for unstable modes (or tachyonic glueballs, in the language of [41]). The point 
of the previous paragraphs is that thermodynamics provides guidance not only on when 
to expect an instability, but also in which mode. 

It is straightforward to start with the lagrangian in (2) and show that linearized 
perturbations to the equations of motion result in the following coupled equations: 

d5F = d*5F + dSif! A *F = 

(22) 

<Vi-16F^F^ = 0. 



Here □ = g^V^dv is the usual scalar laplacian. F in (22) is the background field 
strength in (6): it is the common value of the four Fa- 5F is not the variation in 
F itself; rather, the variation of the Fa is expressed in terms of SF in (21), with 
i = 1. The variation of the field strength is in a direction orthogonal to the background 
field strength of the AdS^TW solution. The graviton decouples from the linearized 
perturbation equations: ST^ V vanishes at linear order in SF because 5F A • Fa = 0. 11 

For comparison, we write down the linearized equations for fluctuations of the other 
scalars: 

□ + A _ 8i^] = (23) 

for i = 2,3. It was shown in [18] that any perturbation involving only matter fields 
satisfying the dominant energy condition could not result in a normalizable unstable 
mode (that is, a normalizable mode which grows exponentially in Lorentzian time). It 
was conjectured [18, 19] that in fact there was no classical instability at all. The scalars 
tfi do not satisfy the dominant energy condition because of the potential term in (2). 
Thus the outcome of our calculations is not fore-ordained by general arguments, and 
we have a truly non-trivial check on the classical stability of highly charged black holes 
in Af = 8 gauged supergravity. In fact, our results turn out to be in conflict with the 
claim of classical stability in [18, 19]. 

11 Besides the all-charges-equal case, we know of one other case where the metric decouples at linear 
order: Qi = Q 3 with Q 2 = Q4 = 0. There may be other cases as well — presumably whenever 
Qa ■ SQa = and 5S = for an unstable eigenvector (SS, SQa) of H^q a . 
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Decoupling the equations in (22) is a chore greatly facilitated by the use of the dyadic 
index formalism introduced in [42]. For the reader interested in the details, we present 
an outline of the derivation in section 4.1. The final result is the fourth order ordinary 
differential equation (ODE) 

\ ' ' (24) 

4Q 2 ( j +9r/9rl^i(r), 

where we have assumed the separated form 5(fi = Ree~ lu)t Ye m S>pi(r), where Y^ m is the 
usual spherical harmonic on S 2 . This is to be compared with the separated equation 
for the other scalars: 



'£ + ft/*-«±il-H + £, rf fc (r) = „ (25) 



for i = 2, 3. 



4.1 Dyadic index derivation of (24) 

To derive (24) using the dyadic index formalism, it is convenient first to switch to 

H signature to avoid sign incompatibilities between the raising and lowering of 

dyadic and vector indices. One introduces a null tetrad of vectors, (/ M , n M , m M , m M ), 
defined so that / M n M = — m M m M = 1 and all other inner products vanish. Next define 



and set D = / M <9 M , A = n M <9 M , 5 = m M <9 M , 5 = m M <9 M . Vector indices are converted into 
dyadic indices by setting v AA = cr^^y^. Dyadic indices are raised and lowered using 
northwest contraction rules with e i = e 01 = Cqi = e 01 = 1. By demanding that a^ A 
is covariantly constant, one can obtain a unique covariant derivative _D M , whose action 
on a spinor is 

D^r = d^ r -ij^ r . (27) 
The so-called spin coefficients, 7 AAsr = a AA^' lSr ' are conven tionally written as 

k e \ ( a (5 



'"^ I < n J Toisr \ J ,,, 
p a\ ( t 7 x 



7l6sr ~U A J 7lisr ~V 7 v 
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A less compressed presentation of dyadic index formalism can be found in [42, 43], and 
the appendix to [44] . 

For AdS^-RN, a convenient choice of the null tetrad and the corresponding nonzero 
spin coefficients are as follows: 



P = (l//, 1,0,0) n" = ^(W, 0,0) 

1 1 

= — =(0,0,1,10300) m M = — ^=(0, 0, 1, -% esc 9) 



(29) 



r 



1 / R cote 

p = — fi= -— 7 = — a = -0= ■=- . 30 

r 2r 4 V8r 



in (29) and (30) we have not yet taken the black brane limit. Taking this limit replaces 
esc by 1 in (29) and sets a = (3 = in (30). Proceeding without the black brane 
limit, we trade the real antisymmetric tensor F^ v for a complex symmetric tensor, 



through the formula 



Ar "U 0) <A 0) ) { ' 



4V2F^ AA a» r = ^Uat + *i°r e Ar ■ (32) 



The factor of 4^/2 in (31) is for convenience: the AdS^KN background has (j)^ = Q/r 2 
and all other components zero. In the same way we trade in 5F^ U for $Ar, whose 
components are 0o, 0i, and 02, with a similar factor of A\/2. Finally, we write (p in 
place of 5(fi to avoid the ambiguity in the meaning of 5. 

The first order equations for the gauge field in (22) can now be cast in dyadic form 
as follows: 



1 

2 

In components, these equations read 



D A r ® Ar + -d AA ^ r e Ar + *i%r) = . (33) 



(D - 2p)0! -(5- 2a)0 o = -<j>^Dip 
(A + n- 2 7 )0 o -^i = O 

(34) 

(D — p)0 2 — 501 = V ' 

(6 + 2/3)0 2 - (A + 2^)0! = 0f } Ay? . 

It is possible to combine these equations into three second order equations in which only 
a single 0, appears. Together with the scalar equation, these equations are equivalent 
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to (22): 



(D - 3p)(A + n - 27) - 5(5 - 2a)] O = -(p^SDp 
"(A + 3n)(D - p) - 5(5 + 2(3)} 2 = -<pf ] 5A v 
(D-2p)(A + 2/i) - {5 + (3-a)5 



□+7> + 2(<M 0) ) 2 



(p = -4^ o) Re0! 



(35) 



where we have made use of the fact that the spin coefficients are all real for AdS 4 - 
RN. The equations (34) are a special case of (3.1)-(3.4) of [45]. The first and second 
equations of (35) are (3.5) and (3.7) of [45], and the third is derived in a similar manner. 
The fourth is the scalar equation in (22), but to preserve the definition of □ we write 

□ = —g^ v \j vi Q v i n _| conventions. The differential operators in the third equation 

of (35) are purely real (this takes a bit of checking for (5 + (5 — a)5), so we can take 
the real and imaginary parts of this equation. The equation for Im0! decouples from 
all the others. The equations for O and 02 are sourced by <p, but 0o and 02 do not 
otherwise enter; thus one can solve first for Re 0i and ip, and afterwards use the first 
and second equations in (35) to obtain O and 02- Since 4>^ is nowhere vanishing, the 
last equation in (35) can be used to eliminate Re 0i algebraically. The final result is 



{D-2p){A + 2p)-{5 + f3-a)5 



1 



>(0) 



□ +— + 2(01' 



(0)\2 



up = <j)f ] DA(p . (36) 



Plugging in the separated ansatz ip = Re {e lu)t Y im 5(pi{r)}, one easily obtains (24). 



4.2 Numerical results from the fourth order equation 

A dynamical instability exists if there is a normalizable, unstable solution to (24) 
or to (25). Neither of these equations admits a solution in closed form, so we have 
resorted to numerics. Briefly, the conclusion is that, in the black brane limit and within 
the limits of numerical accuracy, we find a single unstable mode for (24) precisely 
when x > 1, and no instabilities for (25). This is completely in accord with the 
intuition from thermodynamics: (25) represents a fluctuation that has nothing to do 
with the variation of charges that gave the unstable eigenvector of the Hessian matrix 
of M(S,Qi,Q2,Qs,Qi)- The unstable mode in (24) persists to finite size AdS^-KH 
black holes, but eventually disappears for small enough black holes. 

To carry out a numerical study of (24), the first step is to cast the equation in terms of 
a dimensionless radial variable u, a dimensionless charge parameter x, a dimensionless 
mass parameter a, and a dimensionless frequency Co: 

r Q = (£\" ~ = ^! ( o 7) 

U MV3L2/3 x M 2 / 3 L 1 / 3 ° KM) U MV3 • 
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Then we have 



"~ + d uf d u - 0- 2 ) ( "~ r - °" 2 3 + - T = 

v / u 2 J \ f u 2 u 6 -u 4 y 



^ 2 



4X 2 y + #i (38) 



2 v 2 

/ = (7 - - + — + W 



2 

" T U 

U 



Evidently, the dimensionless control parameters are t (the partial wave number), a, 
and x- Using Mathematica, we solved (38) numerically via a shooting method, and 
obtained wavefunctions 6<pi(r) which fall off like 1/r 2 near the boundary of AdS± and 
at least as fast as (r — rn)^^'^ 11 ^ near the horizon. 

To check that the wavefunction is well behaved near the horizon 12 let us transform 
to Kruskal coordinates. The metric near the horizon is 

dv 2 

ds 2 « -f(r H )(r - r H )dt 2 + //(T , g)(r _ rg) + ^2 , (39) 

where r# is the radius of the horizon. Dropping the S 2 piece and introducing a tortoise 
coordinate r*, null coordinates P±, and Kruskal coordinates (T,R) according to 

dr* 1 



dr f'(r H )(r-r H ) (40) 

p ± = e \f(r H )(±t + r,) = ±T + R ^ 

one finds that the near-horizon metric is indeed regular: 

dr 2 4 

ds\ = -f'(r H )(r - r H )dt 2 + — - — = j^-A-dT 2 + dR 2 ) . (41) 

f'{r H )(r-r H ) f'(r H ) 

Having a radial wavefunction 6<pi(r) = (r — r H ) ( TH ^p(r — r H ), where p(r — r H ) 
remains bounded at the horizon, means that the time- dependent perturbation (with 
angular dependence suppressed) is 

<fy>i(t, r) ~ (r - r H ) H/f ' {rH) e Mt p(r - r H ) ~ P+ M ' // ' (rff) ' p(P + P-) , (42) 

which remains bounded as P_ — > 0. The black hole horizon is at P_ = 0, P + > 
(see figure 3). Thus we see that the perturbation is small at the horizon in good 
coordinates, at least for small P + . (As the perturbation grows, the horizon eventually 
starts to fluctuate, but this is not an issue in the question of whether the instability 
exists). 
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We thank G. Horowitz for suggesting that this check should be made. 
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A qualitative summary of our numerical results is displayed in figure 4(a). An 
example of a normalizable wave-function with negative u 2 is shown in figure 4(b). 
Some points to note are: 

• The boundary of the region of dynamical stability comes from instability in the 
£ = 1 mode. The £ = mode is projected out by charge conservation. Higher £ 
modes become unstable in the upper left part of the shaded triangle in figure 4(b). 
The boundaries of dynamical instability for different £ all come together at a — 0. 

• At a — 0, thermodynamic stability is lost at x — 1, whereas dynamical instability 
sets in at x = 1-007. We believe that the 0.7% discrepancy is due to numerical 
error. 

• We have drawn the regions of dynamical instability and thermodynamic stability 
as disjoint in figure 4(a). In fact, our current numerics shows them overlapping 
by about 0.1% around a = 0.1. We do not view this as significant because the 
numerical errors seem to be around 1%. 

Finally, it is worth pointing out that the string theory program of computing black 
hole entropy via a microscopic state count in a field theory dual (see for example [46], 
or [47] for a review) has proved hard to extend past the boundaries of thermodynamic 
stability. For instance, we have a good understanding of the entropy of near-extremal 
D3-branes [48, 49], but not of small Schwarzschild black holes in AdS. It seems to us 
that this is no accident: most sensible field theories have log-convex partition functions, 
and this translates into Hessian matrices Hf 4 Q A which have no negative eigenvalues. 
Pushing past the boundary of thermodynamic stability in a field theory may be possible 
(particularly as one crosses a phase boundary and begins to nucleate the new phase), 
but doing so seems likely to produce dynamical instabilities in the Lorentzian time- 
evolution. This point of view has indeed informed our entire investigation. 

A dual field theory description of a small Schwarzschild black hole in AdS must in- 
volve thermodynamic instability but no dynamical instabilities. We believe that finite 
volume effects in the field theory are essential in this regard: if one imagines a Landau- 
Ginzburg effective description of the field theory, then derivative terms must restore 
stability to a system whose infrared tendencies are controlled by the thermodynamic 
instability. Various properties of small AdS-Schwarzschild black holes have been ex- 
plored (see for example [50, 51]), but the basic problem of reconciling thermodynamic 
instability with dynamical stability in the presence of a field theory dual remains to be 
addressed. 
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5 Conclusions 



A common conception of the Gregory-Laflamme instability is that a uniform solution 
to Einstein's equations (plus matter) competes with a non-uniform solution, and the 
non-uniform solution sometimes wins out entropically. In such a situation, the generic 
expectation is that there is a first order tunneling transition from the uniform to the 
non-uniform state, which may take place very slowly due to a large energetic barrier. 
In fact, the original papers [1, 2] focused mainly on demonstrating the existence of 
unstable modes in a linearized perturbation analysis of the uniform solution. The 
distinction is between global and local stability. At the level of classical gravity/field 
theory, the latter concept is more meaningful, because with quantum effects suppressed 
it is impossible to tunnel away from a locally stable solution. The aim of this paper and 
its shorter companion paper [17] has been to study local dynamical stability of black 
holes in anti-de Sitter space in relation to a particular notion of local thermodynamic 
stability, namely downward concavity of the entropy as a function of the other extensive 
variables. We reach two main conclusions: 

1. In the limit of large black holes in AdS, dynamical and thermodynamic stability 
coincide. This conclusion is supported by numerical evidence. The small discrep- 
ancy between the observed onset of dynamical and thermodynamic instabilities 
is probably numerical error. 

2. Dynamical instabilities persist for finite size black holes in AdS, down to horizon 
radii on the order of the AdS radius. The evidence is again only numerical, but 
we believe the final answer is correct and robust. 

We regard point 1 as a partial verification of a rather more general conjecture, namely 
that black branes should have Gregory-Laflamme instabilities (in the local, dynamical 
sense of the original papers [1, 2]) precisely when thermodynamic stability is lost. 

Point 2 is surprising because it is the first known example of a stationary black hole 
solution with a point-like singularity which exhibits a dynamical Gregory-Laflamme 
instability. Furthermore, it shows that no-hair theorems cannot always hold in anti-de 
Sitter space. Black branes which experience a Gregory-Laflamme instability are often 
supposed to split their horizons and fall into pieces. For this to happen to a horizon 
whose topology is S 2 and which cloaks a point-like singularity would be truly novel: 
what could tear apart a point-like singularity? 
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Is Cosmic Censorship really threatened by our analysis? 13 It is too early to say. 
Using the heuristic method of calculating the most unstable eigenvector of the Hessian 
of the entropy function, we have argued that adiabatic evolution of unstable black holes 
does not lead to nakedly singular solutions. However this does not bear directly on 
the question of whether the horizon should split as it is assumed to do in the evolution 
of unstable black branes. We leave open many questions as to the eventual fate of 
unstable black holes in AdS: Might they settle down to new non-uniform stationary 
solutions? Do regions of strong curvature form? Does the horizon split? Does AdS 
itself fragment through a classical process? We leave these issues for future work. 
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Figure 1: An example of a mass function whose convex hull is flat. The region we 
interpret as stable is from A to B. 
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(a) (b) 

Figure 2: Plots of the most unstable eigenvector of the Hessian matrix of 
S(M, Qi, Q 2 , Q3, Qa)- The inner curves are boundaries of stability. The outer curves 
(when they are present) denote the boundary between regular black branes and naked 
singularities. 




Figure 3: The Penrose diagram of a regular AdS black hole. We can take T = R = 
P + = P_ = at the center of the diagram. The black hole horizon is the diagonal line 
going up and right from the origin. 
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Figure 4: (a) A topologically correct representation of dynamical and thermodynamic 
stability in the whole \-a plane (but see the text regarding possible overlap of the 
two shaded regions), (b) A sample normalizable wave-function with negative u 2 : here 
a = 0.3, x = 0-96, and Cu 2 = -0.281. 
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